hep-ph/96(JS235 



TECHNION-PH-96-18 
August, 1996 

QCD Corrections to b ^ se^e~ Decay 

Cai-Dian Lii" and Han- Wen Huang^ 
a Department of Physics, Technion-Israel Institute of Technology, 

Haifa 32000, Israel; 
b Institute of Theoretical Physics, Academia Sinica, P.O.Box 2735, 

Beijing 100080, China. 



Abstract 

We give a more complete calculation of 6 ^ se'^e'^ decay including leading log QCD 
corrections from nritap to M-\^ in addition to corrections from Myi/ to rrif,. The differential 
decay rate is found to be slightly suppressed for a large invariant mass of the e'^e~ pair; 
while the integrated width is slightly enhanced comparing with the results without the QCD 
running from ratop to My^. 
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1 Introduction 



The rare decay b —>■ se'^e^ is one of the very useful channels for the study of models beyond 
standard model. This deeply depends on more precise calculations of this decay rate. Although 
some calculations show that there exist large cc resonance contributions, which interfere with 
short distance contributions [0, there is still a window in the invariant mass spectrum of the 
e~^e^ pair for short distance contributions to be dominant Furthermore, window also exists 
for short distance to be dominant in the one lepton energy spectrum 0. The QCD corrected 
coefficients of effective operators from b —>■ se^e~ are also important for the exclusive processes 
such as B ^ K{K*)e+e-. 

The decay of 6 — * se^e~ and its large leading log QCD corrections have already been calculated 
in many papers [^, |^, ^ 0, 0. And also some efforts are made to give a next to leading log 
calculations which is estimated within 20% contribution. All these efforts make it easy for 

experiments to detect this channel. However, all these papers do not include the QCD running 
from rritop to Mw- Since the top quark is found to be 2-times heavier than W gauge boson 
it needs a detail calculation for the effect of the QCD running from m^p to M\y. Further 



more, unlike the b —>■ s'-f case, here the electromagnetic penguin is numerically less important 
than the box diagrams and Z boson penguin, although it is much enhanced by leading-log QCD 
corrections. Meanwhile, the important box diagrams and Z boson penguins have no leading-log 
QCD corrections as a running from M^y to [|12]. The only leading-log QCD corrections may 



come from running from rritop to Mw- 

In the present paper, by using effective field theory formalism in standard model, we recalculate 
the b — > se~^e~ decay including QCD running from rritop to M^, in addition to corrections from 
M\Y to rrib, so as to give a complete leading log results. First in the next section, we integrate 
out the top quark, generating a five-quark effective hamiltonian. By using the renormalization 
group equation, we run the effective field theory down to the W-scale, so as to give out the QCD 
corrections from rritop to M^. In section 3, the weak gauge bosons are integrated out at 
scale. Then we continue running the effective hamiltonian down to b-quark scale to include QCD 
corrections from Mw to m^. In section 4, the differential branching ratio is given as a function 



of the invariant mass of the e"^e pair. Our results will be useful for experiments to distinguish 
backgrounds like cc resonance. Section 5 is a short summary. 



2 QCD Corrections from = rritop to = Mw Scale 

At first, in the standard model Lagrangian, we integrate out the top quark, generating an effective 
theory, introducing dimension-5 and dimension-6 effective operators as to include effects of the 
absent top quark. Higher dimension operators are suppressed by a factor of p'^/mj, where 
characterize the interesting external momentum of b quark ~ m^. For leading order of ml/mf, 
dimension-6 operators are good enough to make a complete basis of operators: 

OIr = -Y^rribSLD^R, 

OIr = l^^'^^^m.-SLa^'^hRF^,, 
Qlr = i^'glmb(l)+(l)-SLbR, 

Pt = ¥'"^,F,-s,Y^'D^h,, 

Rl = il^^gl{D''(f)+)(l)^SLj^bL, 

Wlr = -tl^'glmbWlWt^SLa,,bR, 

Wl = i^x'glWlW^SLl,]PlM. 

Wl = ifi'gliD'^Wl)Wt^SLi,7a7ubL: 

Wl = ifi'glW+^Wi^SLPbL, 

Wt = ilJ^'glWlW^SLiD^l. + lMbL. 

Zl = ^^J''gl -p. — ^-^Z''sL%bL(p'^(f)-, 

COS Bw sm 

Zl = ip^'gl / ■ , Z'^SL^.^^^^bLWi^W^, 
COS U,,, sm Uw 



IGvr^ cos 6,,, sin 6. 



w 



O9 = {eyi6n^){sLYbL)ej^e, 

Oio = (eVl67r2)(sL7'^6L)e7^75e. (1) 
Here 'SLD^'^vbi stands for {siD^^jj^bL + {D^SL)jubL) and the covariant derivative is defined as 



with denoting the QCD couphng constant. The tensor appearing in pI'^ assumes the 



following Lorentz structure, the index A ranging from 1 to 4: 

The subscript L and R in the above formula denote left-handed and right-handed quarks, respec- 
tively. Here operators involving hsz(j)W are not included because their coefficients are suppressed 
by rrib/mt and they do not mix with other operators. Then we can write down our intermediate 
effective hamiltonian: 

Heff = 2V2GFVttVZ J2 (2) 

i 

where Vij represents the 3x3 unitary Kobayashi-Maskawa matrix elements. 

The coefficients Ci{fi = rritop) can be derived through matching Green functions calculated 
from the standard model with that from the intermediate effective theory |jl3[. Keeping only 



leading order oi p^/mf, the coefficients relevant to 6 — S7 decay are already given in ref. |13|, here 
we only give the new ones: 

1 n 2 



1 2 

Czi = — ^-sin^ 9^5, (4) 
93 -J 

1 - 66 1 + 26- 206^ + 126^ 
= m^)' log5-21og5 (5) 

+ sin2^^Q + 25-^log(5-y51og(5), 

1 f 16+^6' S-6' + l6\ \ 

1 [16 + 16' S-6' + l6\ \ 



with 6 = M^r/mf. The renormahzation group equation satisfied by the coefficient functions Ci{fi) 
is 

/^I^C-.l/i) = E(7^)^.C.(/i). (8) 

The anomalous dimension matrix is calculated in practice by requiring renormahzation group 
equations for Green functions with insertions of composite operators to be satisfied order by order 
in perturbation theory. 

Only the last five operators in equation are different from that of 6 — 57 case |]I3|. After 
evaluating the loop diagrams, we get the leading order anomalous dimensions for each of the 
operators in our basis. Since there are so many operators, it is a very large matrix of anomalous 
dimensions [|1^. Here we only list the part new from that of 6 ^ 57 case: 
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487r2 












Zl 
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Zl 




23 
487r2 







Zl 

-| + sin^ Oy, 
—6 + 6 sin^ 6u 
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4 sin^ 8^ 


4 sin^ e 

































(9) 



The solution to renormahzation group equation (H) appears in obvious matrix notation as 



exp / dg—— 



(10) 



After inserting anomalous dimension matrix @, we can have the coefficients of operators at 
II = M^, where the W boson has not been integrated out. 



^ ^ ZZb i 

-^dsi^t) [Cji2^{mt) + 2Cwi{mt)\ log 5, 

Cp.{M+) = Cp.{mt) + 12gl{mt)Cwi{mt)\og6, 
CzAM^) = CCzAmt), ^ = 1,2, 

5 



Cz,iM^) = Cz,{mt) + gl{mt)Cz] {rrit) log 5 + 2gl{mt)Cz^^ {nit) log 5 



Co,(M+) = CoMt) + 



-(^ - sin^ 9^)gl{mt)C^2^{mt) \og6 - 6(1 - sin^ 9^) gl{mt) 0^2 (nit) log 6, 
1 



4 sin^ 9. 



2^glim)Cwf irnt) log 5, 
1 



4r sm f/7/; 



(11) 



with ( = as{mt)/as{Mw)- The coefficients of other operators at yU = M^r like pI'^, pI'^, Wl, R\, 
are given at the appendix of ref . Wi . 



3 QCD Corrections from ii = M^y to fi = Scale 

In order to continue running the operator coefficients down to lower scales, one must integrate out 
the weak gauge bosons W, Z and would-be Goldstone bosons at yU = M]^ scale. The matching 
conditions involving four-quark operators and photon, gluon penguin diagrams are the same as 



6 — s> S7 case[]13[. We here only display the following relations between coefficient functions just 
below(-) and above(+) /i = Mw scale, which is relevant to 6 ^ se^e^: 

Cp^.2{M{y) = Cpi4M+)- 7/9, 

Cp2{M^) = Cp2iM+)~gl{Mw)Cwi{M^) -3/2, 
Cp.{M^) = Cp4(M+) + 9, 

CoJM^^) = Co.„(M-) + ^C.,(Mj,) + j-^ + j-^C(M,,). (12) 

The ffist two terms in the right side of Cog{M-^) and Coio{M-^) arise from box diagrams with W 
boson couples to top quark; the third terms arise from box diagrams with W boson coupling to 
charm quark. The function C{Mw) arise from graphs with a Z gauge boson coupling to the e~^e~ 
pair, while the D{M]v) function arises from graphs with a photon couples to the e+e^ pair. They 
are defined as, 

CiMw) = Czi{M^)-glCziiM+) + i^-sm'9JglCp2^{M+) + 2cos'9^gl^^^^ 



D{Mw) = Cp.{M^)-^Cp.4M^)-^Cp.4M^) + ^Cp.{M^). (13) 
In addition to these, there are new four-quark operators from integrating out the W bosonj^, ||: 

O3 = {sLal^hLa)[{uLpl^iULp) + ■■■ + (hili'^^hLp)], 
Oi = {sLal^hLp)[{uLI3'l^ULa) + ■■■ + ^L^3'l^lhLa)\, 
O5 = {sLal''bLa)[iuRi3'J^URp) + ... + (^Rpl ^^Rp)], 

Oq = {sLc,l^hLp)[{uRp'^^,URa) + ... + (^RplfJ:)Ra)], (14) 

with coefficients 

Ci{Mw)=Q, 2 = 1,3,4,5,6, C2{Mw) = -l. 

There should also be two magnetic moment operators relevant to 6 — > S7, 

^7 = -r^rn^SLa^'^hRF^,, (15) 
Os = -^m,SLa^'''X%RG%, (16) 



their coefficients at /i = M-^ are given in ref. |13 



The operator basis now consists of 10 operators Oi — Oio- The effective hamiltonian appears 
just below the W-scale as 

K// = -j^V,,V,:Y.C,{M^)0,{M^). (17) 

V/ i=i 

If the QCD corrections from m^p to Mw are ignored [by setting asijrit) = as{Mw) in eqn.(0)], 
the above results(|12D would reduce to the previous results 0, ^ exactly, where the top quark and 
W bosons are integrated out together. This is a necessary consistent check. 

The running of the coefficients of operators from fi = Mw to fi = rrib was well described in 
previous papers ^. After running we have the coefficients of operators at ^ = nib scale [1T6 |. 



1=1 



Cio{mb) = Cio(M^), 



asirrib) 



1=1 



(19) 
(20) 



with ?7 = as{Mw)/as{mb), hi,ai and Pi defined in ref.[|T6 



4 Results 



In a spectator model, the inclusive decay B XgC^e' is mainly contributed from the b quark 
decay h se^e^. To leading order, the differential decay rate dT{B Xse'^e~)/ds, where 
"5 = (pt +PeY l^li is given by: 



d 



T{B ^ XcCU) ds 



T{B XsC^e 



a 



QED 



[l + 2s)(\C, 



eff\2 



+ C'- 



+ 4(1 + ^) ICrP + UCyReia 



Here 



r^ff 



Csinib) + g{mjmb, S)(3Ci + C2 + SC-s + C, + SC^ + Cg) 
~g{l, s)(4C3 + 4C4 + 3^5 + Ce) - ^g{0, s){C3 + 3C4), 



(21) 

(22) 
(23) 



where g{mc/mfy, s) arises from the one-loop matrix element of the four- quark operators, which 
can be written as 

^1-/3+1 

8 , 8 4^2, n I " 



'1-13-1 



/3<1, 



(24) 



2arctan(l/v/;5^), (3 > I, 
with P = iz^/s. The factor of f{mc/mh) arises because of the dependence of the semileptonic 
decay rate on the nonnegligible ratio mc/rrih, and is given by 

f(x) = l-8x'^ + 8x^ -x^ - 24x^ In x. 

In eqn.(^) the masses of the electron and the strange quark were neglected. 

The semileptonic decay B XceV is used to eliminate large uncertainties of in the 
decay width formula. The dependence on the weak mixing angles also cancels out. If we take 
experimental result Br{B XcCV) = 10.8% [1^, the branching ratios of -B ^ Xse~^e~ is found. 
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Taking values as Mw = 80.22GeV, rrit = 175GeV, rrib = 4.8GeV and the QCD coupling 
constant as{mz) = 0.117 [|I5|, the differential branching ratios are depicted in Fig.l as a function 
of s. It is shown that the QCD corrections from rrit to Mw slightly suppress the b se^e~ 
differential decay rate for a large invariant mass of the e^e~ pair; while enhance it for a very 
small invariant mass of the e~^e~ pair. Since the enhancement is larger than the suppression in 
the invariant mass spectrum of the e+e^ pair, the total branching ratio is found to be slightly 
enhanced about 4% comparing with the one without QCD correction from rritop to M]y 

5 Conclusion 

As a conclusion, we have given the full leading log QCD corrections (including QCD running from 
rritop to Mw) to 6 — se^e~ decay in the standard model. 

The QCD correction from rrit to Mw slightly suppresses the b se^e~ differential decay 
rate for a large invariant mass of the e+e" pair. While the integrated width is slightly enhanced 
comparing with that without the QCD running from rritop to Mw- 

Although this result is not quite different from the previous calculations, our improvement lies 
in reducing some theoretical uncertainties. 
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Figure Captions 

Fig.l Differential branching ratios of 6 — > function of s = {p^ +p Y /ml. The sohd hne 

denotes the result with full QCD corrections, while the dashed one corresponds to result without 
QCD corrections from rrit to Mw- 
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